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青藏 高 原 地 区 山 生 柳 遗 传 多 样 性 研究 


R, L EAR! AKRE 23， 张 发 起 1， 陈 世 龙 UU 

(1. 中 国 科学 院 高 原生 物 适 应 与 进化 重点 实验 室 ， 中 国 科 学 院 西北 高 原生 物 研究 所 ， 西 5 
810001; 2. 中 国 科学 院 大 学 , 北京 100039; 3. 青海 省 作物 分 子 育种 重点 实验 室 , 西宁 810001) 
摘要 ， 利用 叶绿体 非 编码 区 片段 研究 分 布 于 青藏 高 原 地 区 的 山 生 柳 居 群 遗 传 多 样 性 ， 对 未 来 山 
生 柳 生态 环境 的 保护 以 及 对 青藏 高 原 地 区 物种 丰富 度 的 保护 都 具有 指导 意义 。 设 计 并 筛选 出 
cpDNA 引物 5'trnG2G ~ 3'trnG CUUC) 和 5^rpSI2 ~rpL20， 用 扩 增 出 的 片段 和 对 应 的 联合 片段 
进行 后 续 的 遗传 多 样 性 分 析 。 山 生 柳 的 联合 片段 检测 到 35 种 单 倍 型 ， 单 倍 型 多 态 性 0.626， 核 
苷 酸 多 态 性 0.00085。 中 性 检验 Tajima's D (-2.287, P « 0.01) 和 Fu's Fs (-5.295, P < 0.02) 都 是 
著 负 值 ， 推 测 山 生 柳 个 体 数 近期 经 历 过 扩张 。AMOVA 分 析 显 示 ， 居 群 内 和 居 群 间 遗 传 变异 
对 应 93.70% 和 6.30%， 表 明 居 群 内 的 变异 是 山 生 柳 遗 传 变异 的 主要 来 源 。 居 群 间 遗传 分 化 程度 
中 等 偏 低 (Fsr = 0.063)， 基 因 流 (Nm) 为 7.439， 说 明山 生 柳 各 居 群 的 基因 交流 非常 频繁 ， 
同 地 理 居 群 间 存 在 一 定 的 基因 流动 。 遗 传 分 化 系数 Nsr (0.075) 大 于 Gsr (0.068). 和 基于 遗传 
型 的 UPGMA 聚 类 分 析 , 表明 山 生 柳 12 个 居 群 分 为 4 组 且 与 居 群 的 地 理 分 布 没有 明 
相关 性 。 山 生 柳 是 进行 有 性 繁殖 还 是 无 性 繁殖 主要 受 环境 因素 的 影响 ， 居 群 内 变异 是 山 生 柳 
遗传 变异 性 的 主要 来 源 ， 居 群 间 基因 交流 频繁 。 
关键 词 : 山 生 柳 ， 基 因 流 ， 青 藏 高 原 ， 叶 绿 体 DNA， 遗 传 多 样 性 
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Genetic diversity of Salix oritrepha (Salicaceae) in the 
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Abstract: Using the non-coding region of chloroplast DNA fragments to study the genetic diversity of 
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Salix oritrepha in the Qinghai-Tibetan Plateau. The study has guiding significance for the protection of 
ecological environment and animal husbandry of S. oritrepha and species richness in Qinghai-Tibet 
Plateau in the future. We designed two useful primers 5'trnG2G ~ 3'trnG (UUC) and 5'rpS12 ~ rpL20. 
The amplified fragments and their combined fragments were used for subsequent genetic diversity 
analysis of S. oritrepha. We detected 35 haplotypes by the combined sequences, haplotype diversity 
(Hd) is 0.626, nucleotide diversity (Pi) is 0.00085. Neutrality tests of Tajima's D (-2.28670, P « 0.01) 
and Fu's Fs (-5.29805, P « 0.02) showed significantly negative values, it is means that S. oritrepha 
had expansion recently of its individual numbers. Analysis of molecular variance showed that 
within-population variance was described as 93.70% of the total variance, while among-population 
variance is 6.3095, it shows that within-population variance is the main resource of genetic variation 
of S. oritrepha. The fixation index (sr) is 0.063 and number of migrants (Nm) is 7.439, it means that 
gene flow is very frequent among different geographical populations. An estimation of 
non-signification Nsr (0.075) is bigger than gene differentiation coefficient CsT (0.068). UPGMA tree 
based on genetic distance and haplotypes suggested that all populations were divided into 4 groups and 
had no significant correlation with the geographical distribution of the population. The mode of 
reproduction (sexual or asexual) of S. oritrepha is mainly depends on environmental factors. The 
genetic variability mainly comes from within- population and it has a frequently gene flow among 
populations. 


Key words: Salix oritrepha, gene flow, Qinghai-Tibetan Plateau, chloroplast DNA, genetic diversity 


ÆI (Salix oritrepha) 是 中 国 特有 物种 ， 落 叶 矮 小 灌木 ， 高 10 ~ 120 cm， 藏 药 中 称 朗 玛 ， 
生 于 海拔 3000 ~ 4700m 的 高 山地 带 ， 繁 衍 方式 以 有 性 生殖 为 主 ， 无 性 繁殖 为 辅 ， 是 青藏 高 原 地 
区 高 寒 柳 灌 丛 的 建 群 种 和 优势 种 ， 具 有 极 高 的 药 用 价值 ， 对 生态 环境 的 保护 也 尤为 主要 。 一 方 
面 对 山 生 柳 化 学 成 分 的 研究 ， 证 明 其 可 解 热 、 消 炎 止 痛 、 治 疗 妇科 疾病 和 提高 机 体 免疫 力 〈 封 
士 兰 等 ，2001)。 男 一 方面 山 生 柳 是 深 根 性 灌木 ， 根 深长 度 达 到 植株 纵向 高 度 的 4~5 倍 ， 故 可 
以 抵抗 干旱 和 低温 。 由 于 耐 高 盐 高 碱 且 生长 繁殖 快速 ， 故 能 涵养 水 源 、 治 沙 和 防 沙 ， 使 其 起 到 
河岸 防护 、 水 土 保持 的 作用 ( 陈 文 业 等 ，2008)。 常 与 高 山 绣 线 菊 (Spiraea alpina)、 金 露 梅 
(Potentilla fruticosa), Vb3& (Hippophae rhamnoides), U EIE CSibiraea angustata), KE 
WI (Salix pseudospissa) 等 组 成 高 寒 江 从“〈 王 芳 等 ，2010)， 是 藏 民 牲 畜 的 夏季 放牧 地 。 

遗传 多 样 性 分 析 是 评价 和 保护 濒危 物种 的 主要 指标 〈Cires et al., 2011)， 种 内 遗传 多 样 性 若 
愈 丰 富 ， 物 种 对 环境 的 适应 能 力 就 愈 强 ,保护 了 物种 和 生态 系统 的 多 样 化 ， 也 能 减 慢 因为 适应 、 
进化 所 致 的 灭亡 进程 。 近 几 十 年 来 , 对 青藏 高 原 地 区 植被 的 遗传 多 样 性 的 研究 越 来 越 深 入 (Duan 
et al., 2011; Khan et al., 2014)， 这 些 都 对 保护 青藏 高 原 地 区 植被 丰富 度 具 有 重要 的 指导 意义 。 研 
究 植 物 遗 传 多 样 性 的 方法 中 ,叶绿体 基因 组 DNA. 被 广泛 使 用 (Demesure et al., 1996; Rungis et al., 
2017; Wu, 2016)。 叶 绿 体 DNA 的 结构 紧凑 简单 、 分 子 量 小 ， 可 从 分 子 角度 作 比较 精准 的 分 析 ， 
其 比 核 DNA 存在 更 加 稳定 、 进 化 更 加 缓慢 且 编 码 序列 也 保守 ， 故 广泛 应 用 于 植物 遗传 多 样 性 的 
剖析 研究 (Shaw et al., 2014)。 叶 绿 体 DNA 在 大 多 数 被 子 植物 中 为 母系 遗传 ， 只 通过 种 子 进行 
遗传 ， 不 受 花粉 遗传 的 影响 (Mogensen, 1996). 
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目前 对 高 寒山 生 柳 灌 从 的 群落 特征 的 研究 、 在 畜牧 业 中 的 作用 以 及 在 防止 沙漠 化 中 起 到 的 
作用 都 有 少量 研究 (Kang et al., 2015; 王 纬 和 王 盛 ，2015; 杨 路 存 ，2013)， 但 对 山 生 柳 遗 传 多 
样 性 的 研究 较 少 ， 之 前 的 研究 表明 山 生 柳 种 群 具 有 较 高 的 遗传 多 样 性 且 随 着 海拔 梯度 的 升 高 而 
呈现 升 高 趋势 〈 郭 敏 等 ，2012)， 并 且 遗 传 变 异 主要 存在 于 居 群 内 。 由 于 其 所 研究 的 地 理 分 布 范 
围 仅 限 祁连山 ， 不 能 准确 分 析出 山 生 柳 遗 传 结构 及 种 群 动态 〈 郭 敏 等 ，2012 )。 

本 研究 采集 了 青藏 高 原 大 部 分 地 区 的 山 生 柳 ， 设 计 并 筛选 出 两 对 特异 性 针对 山 生 柳 的 
cpDNA 引物 ， 用 叶绿体 非 编码 区 的 联合 片段 对 山 生 柳 进 行 遗传 多 样 性 的 研究 。 本 研究 将 对 山 生 
柳生 殖 系统 的 研究 、 种 质 资源 的 保护 和 生态 环境 的 保护 都 有 重要 意义 。 


1 材料 与 方法 


1.1 材料 


2006 年 ~ 2017 年 对 分 布 于 青藏 高 原 地 区 的 山 生 柳 进行 广泛 的 野外 采集 。 对 每 个 野外 自然 
X 居 群 的 选取 8 ~ 58 个 个 体 ， 个 体 间 距 超 过 100 m， 共 收集 山 生 柳 12 个 居 群 351 个 个 体 〈 表 D. 
© 采集 山 生 柳生 长 良好 的 新 鲜 的 幼 嫩 叶 片 ， 迅 速 用 硅胶 来 进行 干燥 ， 而 后 放 于 -20'C 长 期 保存 。 所 
e 采集 的 凭证 标本 存放 于 中 国 科学 院 西北 高 原生 物 研究 所 青藏 高 原生 物 标本 馆 (HNWP)。 

表 1 山 生 柳 各 居 群 采集 信息 


Table 1 The information of 12 populations of S. oritrepha 


编号 凭证 标本 采集 地 纬度 NO 经 度 (E) 海拔 (m) ”个 体 数 
Code Voucher specimen Locality Latitude Longitude Altitude Number 
GZ Chen2006247 四 川 ， 甘 孜 Ganzi, SC 29"00' 99°44" 4170 27 
MQ Chen2006328 青海 ， 玛 沁 Maqi, QH 3418' 10017' 3990 29 
YS zhang2016713 青海 ， 玉 树 Yushu, QH 32/53! 9634' 4425 57 
g QM chen2014621 海 ， 曲 麻 莱 Qumalai, QH 33°12' 9727 4440 8 
GD zhang2016125 青海 ， 贵 德 Guide, QH 3618' 10138' 3103 43 
MY Gao2017087 青海 ， 门 源 Menyuan, QH 37 48' 101708' 3520 58 
DR chen2013173 青海 ， 达 日 Dari, OH 3318' 100727' 4130 15 
BM zhang2015118 青海 ， 班 玛 Banma, QH 32°42' 100°40' 4464 28 
NQ zhang2016468 THÉ. SÉWÉNangqi, QH 31'00' 9622! 4085 34 
QL Gao2017016 青海 ， 祁 连 Qilian, QH 38'45' 98^12' 3786 24 
SN Gao2017044 甘肃 ， 肃 南 Shunan, GS 3843' 99728' 2901 13 
DT Gao2017104 ij. A38 Datong, QH 37T12' 101 48' 3681 15 


1.2 DNA 提取 与 PCR 扩 增 


选取 改良 的 CTAB 法 从 干燥 处 理 的 山 生 柳 嫩 叶 中 提取 基因 组 DNA CPorebski et al., 1997). 
参照 序列 选取 已 上 传 到 NCBI 上 的 Salix arbutifolia, Salix babylonica , Salix interior ~ Salix 
oreinoma . Salix purpurea. Salix suchowensis 和 Salix tetrasperma 等 柳 属 植物 叶绿体 基因 组 序列 ， 
引物 片段 的 选取 参考 已 发 表 过 的 期 刊 文章 (Hamilton, 1999; Scarcelli et al., 2011; Shaw et al., 


2005), 


J A 


F Primer Premier 5.0 (Singh et al., 1998), BioXM 2.6 CRIE PIKT 4 


> 2004) 和 


MT 


SnapGene 3.2.1 (McKeone et al., 2014) 进行 引物 设计 ， 共 设计 了 21 对 引物 用 于 PCR $735. 37 
增 反 应 体系 为 25 uL: 2.5 uL 的 10xPCR Buffer C 1.5 mmol-L! MgCb), 0.5 uL 的 10 mmol LE! 
dNTP, 1E 5| 7] 0.5 uL C10 pmol-L-D, 1 个 单位 的 Taq DNA 酶 (TaKaRa, KÆ), 15 — 50 ng 
的 模板 DNA， 去 离子 水 补 齐 到 25 uL. PCR 反应 扩 增 程序 : 94 C 预 变性 3 min; 94 C 变性 45 
s; 60 COB/AX 45s; 72 CÁE[f 30s; 35 个 循环 ， 72 C Ef 5 min; 4 C 保存 。PCR 反应 产物 
送 至 生物 工程 (上 海 ) 股份 有 限 公 司 测序 。 将 获得 21 对 引物 在 12 个 居 群 24 个 个 体 中 进行 PCR 
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获得 两 对 遗传 多 态 性 较 高 的 引物 [5’trnG2G ~ 3'tmG. (UUC) 和 S"rpS12 ~ 


rpL20] 用 于 所 有 个 体 的 扩 增 测序 〈 表 2)。 
表 2 用 于 PCR 扩 增 的 cpDNA 引物 序列 
Table 2 The cpDNA primer sequences for PCR amplification 


编号 引物 名 称 退火 温度 引物 序列 (5'~ 3) 
Code Primer name Temperature Primer sequence (5'~ 3') 
1 5'trnG2G-3'trnG (UUC) 60 °C F CGAACCCGCATCATTAGCTT 

R GCAGCATAGCGGGTATAGTT 
2 S’rpS12-rpL20 60 °C F ACCAAATCCCCCGCTCT 


R GCCCCACGAAAGCTTGATGT 


1.3 数据 分 析 


运用 Chromas 软 伯 
2018) 对 所 有 测序 序列 进行 人 了 


位 排列 。 利 
Tajima's D 和 


方法 (analysis of molecular variance, AMOVAO ) 检测 居 群 内 和 居 群 间 的 遗传 变异 水 平 
型 分 布 的 Fsr 进 行 评价 〈1000 2X8 


关系 计算 Nm 值 


并 用 U- 统 计 方法 对 Gsr 和 Nsr 进行 比较 (1000 次 重复 置换 检测 ) 以 检测 山 生 柳 单 1 
点 (Caraux & Pinloche, 2005 )。 利 
TCS 简约 网 络 关 系 图 。 


F Chttp://technelysium.com.au/wp/chromas/) 和 MEGA 7.0.26 (Kumar et al., 
[校对 ， 并 用 MEGA 软件 对 校对 好 的 所 有 个 体 的 扩 增 片段 进行 对 
] DnaSP 6.10.04 (Rozas et al., 2017). 软件 确定 山 生 柳 单 倍 型 ， 并 统计 中 性 检验 中 的 
Fu's Fs。 运 用 SPSS (George & Mallery, 2003) 软件 获得 基于 遗传 距离 和 单 
布 的 UPGMA 聚 类 


倍 型 分 


分 析 图 。 利 用 Arlequin 软件 包 CExcoffier & Lischer, 2010) 中 的 分 子 变异 分 析 


并 对 单 售 


EE 复 置换 检测 )。 根 据 基因 流 Nm) 与 遗传 分 化 系数 Fsr 之 间 的 


， 居 群 平均 等 位 基因 数 。 失 配 分 布 (mismatch distribution?) 分 析 也 在 Arlequin 
中 完成 。 通 过 PERMUT 软件 计算 居 群 内 平均 遗传 多 样 性 
populations, As) ~ ME Z PETE C 


(average gene diversity within 


居 群 件 遗 传 分 化 系数 (Gsr) 和 (Nsr)， 


total gene diversity, hr) ~ 


型 的 分 布 特 


OK 


2 结果 与 分 析 


] TCS 软件 (Clement et al., 2000) 来 建立 山 生 柳 单 倍 型 之 间 的 


2.4 山 生 柳 遗传 多 样 性 分 析 

对 351 个 个 体 的 扩 增 产物 进行 单 向 测序 ， 用 MEGA 软件 对 位 排列 且 删 除 poly 结构 后 ， 
5'trnG2G ~ 3'trnG. (UUC) 和 5’rpS12 ~ rpL20 片段 矩阵 长 度 分 别 为 481bp、730bp， 共 检测 到 36 
个 变异 位 点 ， 其 中 17 个 为 简约 性 信息 位 点 ， 变 异 率 为 2.97%。 山 生 柳 的 联合 片段 检测 出 35 个 
单 倍 型 (Haplotype, H01 ~ H35), 5'trnG2G ~ 3'trnG (UUC) 片段 的 单 倍 型 序列 和 5’rpS12 ~ rpL20 


片段 的 单 倍 型 序列 提交 至 GeneBank (MK602662 ~ MK602696 和 MK606373 ~ MK606407)。 联 
合 片 段 的 单 倍 型 多 样 性 0.6261， 单 倍 型 多 样 性 标准 差 0.029( 表 3)。 各 居 群 联合 片段 之 间 的 多 态 
性 表现 不 同 ， 居 群 的 平均 等 位 基因 数 (Na) 依次 为 1.067、1.078、1.278、1.022、1.656、1.089、 
1.067、1.100、1.056、1.078、1.089、1.033、1.134。 遗 传 多 样 性 最 丰富 的 是 居 群 YS， 其 具有 最 
高 的 核 昔 酸 多 样 性 〈0.00145)、 平 均 核 音 酸 差异 〈1.722) 和 最 多 的 单 倍 型 数目 (132 CE 4D. 
表 3 山 生 柳 单 倍 型 多 样 性 分 析 和 中 性 检验 

Table 3 The haplotype diversity analysis and neutrality test of of S. oritrepha 


基因 片段 单 倍 型 单 倍 型 单 倍 型 多 样 性 单 倍 型 多 样 性 中 性 检验 中 性 检验 
多 样 性 (Hd) 数目 (h) HÆ Wh) 标准 差 (Sh) 
Fragment Haplotype No. of Variance of Standard deviation Tajima’s D Fu's Fs 
diversity Haplotype haplotype of haplotype 
diversity diversity 
5’trnG2G ~ 3'trnG -2.46738 -7.06562 
0.2291 22 0.0008 0.02 
(UUC) i 2 P «0.001 P< 0.02 
-1.70183 -2.22048 
S'rpS12 ~ rpL20 0.5472 18 0.000 0.02 
人 "n s 0.05-P«0.10  0.05<P<0.10 
联合 序列 -2.28670 -5.29805 
0.6261 35 0.00085 0.029 
Combined segment P «0.01 P«0.02 


由 各 居 群 的 单 倍 型 分 布 可 知 ， 仅 H01 被 全 部 12 个 居 群 所 共享 ， HO4. H05, H06 和 H07 是 
GZ 所 拥有 的 特有 单 倍 型 ，H11 是 MQ 所 拥有 的 特有 单 倍 型 ，H17、H18、H19、H20、H21 和 
H32 是 YS 所 拥有 的 特有 单 倍 型 ，H27、H28、H29 和 H30 是 GD 所 拥有 的 特有 单 倍 型 ，H16 是 
MY 所 拥有 的 特有 单 倍 型 ，H15 是 DR 所 拥有 的 特有 单 倍 型 ，H24、H25 和 H26 是 BM 所 拥有 的 
特有 单 倍 型 ，H31 是 NQ 所 拥有 的 特有 单 倍 型 ，H33 是 QL 所 拥有 的 特有 单 倍 型 ，H34 和 H35 
是 SN 所 拥有 的 特有 单 倍 型 ， H02，H03，H08，H09，H10，H12，H13，H14，H16，H17，H22， 
H23 是 2~ 10 个 居 群 所 共有 的 单 倍 型 《图 1 和 表 4) 
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山 生 柳 居 和 群 和 单 倍 型 的 分 布 图 


Distribution of populations and haplotypes of S. oritrepha 


倍 型 Haplotype 的 简写 。 
Note: GZ, MQ, YS, OM, GD, MY, DR, BM, NQ, QL, SN, DT are the code name of 12 populations; 
“Hap” is the abbreviation of Haplotype. 


E: GZ, MQ, YS, QM, GD, MY, DR, BM, NQ, QL, SN, DT 均 为 12 居 群 的 编码 名 称 ; 
“Hap” 73 ff. 


X A 山 生 柳 居 群 单 倍 型 多 样 性 分 布 


Table 4 The haplotype diversity analysis of all populations of S. oritrepha 


居 群 单 倍 型 多 样 性 CH4) PARAH OD 单 倍 型 分 布 
Population Haplotype diversity No. of Haplotype Distribution of Haplotypes 

GZ 0.593 7 H01 (17), H02 (04), H03 (01), H04 (02), H05 (01), 
H06 (01), H07 (01) 

MQ 0.798 7 H01 (11), H02 (06), H08 (01), H09 (04), H10 (03), 
H11 (03), H12 (01) 

YS 0.626 13 H01 (32), H02 (09), H03 (01), H09 (03), H13 (01), 
H14 (01), H17 (02), H18 (01), H19 (01), H20 (01), 
H21 (01), H22 (02), H32 (02) 

QM 0.464 3 H01 (06), H22 (01), H23 (01) 

GD 0.521 10 H01 (29), H02 (03), H10 (01), H12 (01), H13 (03), 
H23 (01), H27 (02), H28 (01), H29 (01), H30 (01) 

MY 0.589 8 H01 (38), H02 (04), H03 (01), H09 (04), H10 (02), 
H12 (02), H13 (05), H16 (02) 

DR 0.562 5 H01 (10), H02 (01), H09 (02), H14 (01), H15 (01) 

BM 0.265 5 H01 (21), H02 (01), H24 (03), H25 (02), H26 (01) 


NQ 


QL 


SN 


DT 


0.701 


0.622 


0.667 


0.371 


HO1 (11), H02 (15), H03 (01), H09 (04), H10 (01), 
H31 (02) 
H01 (14), H02 (03), H08 (01), H09 (02), H12 (01), 
H13 (02), H33 (01) 
HO1 (07), H09 (02), H10 (01), H12 (01), H34 CODD, 
H35 (01) 
H01 (12), H02 (01), H09 (01), H23 (01) 


2.2 山 生 柳 遗传 结构 分 析 


分 子 变 异 分 析 (AMOVA ) 结果 表明 , 居 群 间 遗 传 变异 为 6.30%, 
(P<0.01)， 说 明 居 群 内 遗传 变异 为 山 
居 群 间 遗 传 分 化 很 小 。 根据 基 


0.06298, 


生 柳 遗传 变异 的 3 
因 流 Nm 与 Fsr 之 | 


居 群 内 遗传 变异 为 93.70%， 
EHK CR $)。 遗 传 分 化 系数 Fs 
目的 关系 《Kirk & Freeland, 20112, Nm 


= (1-Fsr) /2Fsr， 得 到 Nm 为 7.439， 基 因 流 很 强 。 软 件 PERMUT 和 DnaSP 算得 ， 山 生 柳 居 群 


内 平均 遗传 多 样 性 (Hs) 为 0.598 (0.0370)、 总 的 遗传 多 样 改 


遗传 分 化 Gsr=0.068，Nsr 70.075 CNsr» Gsr,. P» 0.057. 


E (Hr) 23 0.646 (0.0401), JE HE 


X5 基于 联合 片段 的 山 生 柳 分 子 变异 分 析 


Ei 


Table 5 The AMOVA analysis based on combined sequences of S. oritrepha 


变异 来 源 自由 度 (df) 总 方差 (SS) 变异 成 分 变异 比例 C90) 
Source of variation Degrees of freedom Sum of squares Variation component Percentage of variation 
居 群 间 
11 16.697 0.02586 Va 6.30 
Among population 
居 群 内 
339 264.905 0.78143 Vb 93.70 
Within population 
总 计 Total 350 281.601 0.80729 100 
固定 系数 Fixation index (Fsr) Fsr = 0.06298 


山 生 柳 各 居 群 间 的 遗传 距离 大 部 分 都 较 小 , 其 中 QM 和 NQ 之 间 昌 
YS 和 SN 出 现 了 最 小 遗传 距离 0.000， 大 多 数 
距离 为 0.424; 居 群 间 的 遗传 分 化 指数 CFsr) 显 


间 出 现 了 最 大 的 Fsr = 0.2550 CR 6)。 


bb 现 了 最 大 遗传 距离 0.977， 


示 有 18 个 负 值 ， 


居 群 之 间 的 遗传 距离 差异 性 不 大 ， 总 的 平均 遗传 
H 27.28%。 其 中 QM 和 NQ 之 


表 6 山 生 柳 各 居 群 遗传 距离 CE) 和 分 化 指数 (下 ) 
Table 6 The genetic distance and the fixation indexes of different populations of S. oritrepha 

GZ MQ YS QM GD MY DR BM NQ QL SN DT 
GZ 0.431 0.129 | 0.79 | 0.5004 0244 0315 . 0252 0.550 0.233 0.191 0.235 
MQ 0.0249 0.263 0.765 | 0438 . 0261 0431 0.27 0.534 0.289 — 0.228 — 0331 
YS 0.0057 0.0230 0.682 — 0.2285 0.059 0.106 — 0.072 0.354 0.073 | 0.000 — 0.044 
QM 0.0910 0.1085  -0.0088 0.906 — 0.617 0.881 . 0.676 0.977 0.888 0.704 . 0.673 
GD 0.0089 0.0261 0.0148 -0.0394 0.284 0.473 0.405 0.557 0.316 | 0279 . 0419 
MY 0.0135 0.0304 0.0035 0.0432 0.0150 0.238 . 0.134 0.375 0.120 — 0.078 — 0.169 
DR 0.0882 0.0238 -0.0122 0.0034 -0.0134 -0.0091 0.240 0.548 0.250 — 0.149 — 0.163 
BM 0.0089 0.1229 0.0413 0.0044 0.0262 0.0849 0.0385 0.407 0.126 — 0.108 — 0.185 
NQ 0.0624 0.0074 0.0695 0.2550 0.0590 0.0027 0.1289 02221 0.396 — 0.358 | 0.397 
QL 0.0123 0.0310 -0.0095 0.0245 -0.0041  -0.0191 -0.0206 0.0582 0.1024 0.102 — 0.184 


SN 0.0619 0.0446 0.0040 -0.0086 -0.0112 0.0166 -0.0523 0.0370 ^ 0.1585 -0.0011 0.075 
DT 0.0669 0.0891 -0.0080 -0.0324 -0.0156 0.0133 -0.0237 0.0252 0.2154 -0.0040 -0.0154 
注 : 遗传 分 化 指数 (Fsr)〔 对 角 线 下 方 );， 遗传 距离 (对 角 线 上 方 )。 


Note: Fixation indexes of genetic variation (below diagonal); 


山 生 柳 12 个 居 群 的 遗 


传 距离 的 范 


围 为 0.000 ~ 1.234， 表 明山 生 柳 


the genetic distance (above diagonal). 


导 群 内 遗传 多 样 性 丰富 。 


基于 遗传 距离 的 UPGMA 聚 类 分 析 , 得 到 山 生 柳 的 12 个 居 群 分 为 I(YS, SN, MY, DT, QL, BM), 
II (GZ, MQ, DR), III (GD, NQ) 和 IV (QM) 共 四 组 (图 2)。 基 于 单 倍 型 的 UPGMA 聚 类 分 


H, 得 到 12 个 居 群 分 为 1 CYS, GD, MY), II (NQ), III (BM, GZ) £l IV COM, SN, OM, QL, DT, 
DR) 共 四 组 (图 3)。 两 者 均 没 有 表现 出 与 12 个 居 群 的 地 理 分 布 相对 应 的 聚 类 关系 。 
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Fig.2  UPGMA tree based on the genetic distances 


es | 


遗传 距离 的 UPGMA 


HEZ H 
聚 类 


3 于 单 倍 型 的 UPGMA 


Fig.3 UPGMA tree based on haplotypes 


ES 


2.3 山 生 柳 种 群 动态 
中 性 检验 的 结果 表明 ， 山 生 柳 各 个 引物 片段 和 联合 片段 的 Tajima's D (-2.28670, P < 0.01) 
和 Fu's Fs (-5.29805, P < 0.02) 均 为 显著 的 负 值 〈 表 3)。 失 配 分 析 显 示 ， 在 山 生 柳 分 布 范围 内 


其 分 布 曲 线 整体 呈现 单 峰 分 布 ， 其 SSD Csum of square deviations) f& (0.0038, P — 0.68) 和 Hrag 


N 


(Harpending's raggedness index) 值 (0.0479, P — 0.720 均 不 显著 (P> 0.05)， 这 与 中 性 检验 的 


结果 相 一 臻 《图 4)。 


E E Freq. Obs. 


0.3 
一 一 = HIS Freq. Exp. 


Sis Frequency 


1 2 3 4 5 6 J 8 9 10 11 12 13 14 15 16 17 
失 配 分 析 Pai rwise mismatch 


4 基于 山 生 柳 单 倍 型 的 失 配 分 析 


Fig.4 Mismatch distribution of haplotypes of S. oritrepha 
用 山 生 柳 叶绿体 基因 所 检测 出 的 35 个 单 倍 型 ， 建 立 单 倍 型 多 样 性 的 简约 网 络 图 〈 图 5)。 回 


L 


Hz EHE R, XEJETEINAIH AG CURRO 单 倍 型 即 HO CSaltonstall, 20025. d 
略 去 了 单 倍 型 的 频率 数值 。 图 中 “o” 表 示 理 论 分 析出 的 、 但 没有 实际 观测 到 的 单 倍 型 。 除 YS 
(H32) 和 DG (H27, H29) 的 单 倍 型 同 其 它 单 倍 型 关系 较 远 外 ， 大 部 分 单 倍 型 之 间 仅 有 一 到 
两 步 突 变 ， 演 化 关系 呈现 星 状 辐射 结构 〈 杨 玉 慧 等 ，2004)。 
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H29 H32 
o— 
图 5 山 生 柳 叶 绿 体 基因 单 倍 型 多 样 性 简约 网 络 图 
Fig.5 Parsimony network of chloroplast haplotype diversity of S. oritrepha 
3 讨论 


3.1 山 生 柳 的 居 群 遗传 多 样 性 


筛选 出 的 山 生 柳 特 异性 条 带 为 山 生 柳 无 性 系 鉴定 提供 可 操作 性 。 山 生 柳 的 联合 片段 共 检 测 
到 36 个 变异 位 点 和 17 个 简约 性 的 信息 位 点 ， 插 入 缺失 变异 位 点 19 个 ， 占 总 变异 的 52.7896. 
插入 缺失 是 山 生 柳 变异 位 点 出 现 的 主要 来 源 ， 这 与 DNA 在 进行 复制 时 存在 高 频率 的 点 突变 有 
关 。 在 残酷 的 高 山 山地 冰雪 寒 温 带 生 存 环境 中 , 山 生 柳 的 DNA 复制 发 生 碱 基 错 配 的 可 能 性 更 大 ， 
机 体 的 自我 修复 功能 会 受到 一 定 阻碍 , 说 明 恶 劣 的 生存 环境 可 能 在 一 定 程度 上 会 影响 DNA 复制 
过 程 (Baccarelli & Bollati, 2009; Chinnusamy & Zhu, 2009; Dubrova, 2003). 

对 山 生 柳 居 群 的 联合 片段 进行 数据 分 析 ， 检 测 到 35 个 单 倍 型 。 数 据 分 析 的 结果 总 体 呈 现 出 
较 低 的 遗传 多 样 性 (Hd = 0.6261, Pi = 0.00085)。 山 生 柳 每 个 居 群 间 的 单 倍 型 多 样 性 (Hd) 表现 
不 同 ， 推 测 可 能 是 由 于 在 青藏 高 原 地 区 的 山 生 柳 分 布 范围 广 ， 各 居 群 所 在 的 自然 生存 环境 条 具 
有 差异 ， 经 过 长 期 不 同 的 生态 环境 和 自然 选择 ， 会 导致 具有 一 定 地 理 隔离 的 山 生 柳 居 群 存在 差 
异 的 遗传 多 样 性 。 根据 之 前 对 山 生 柳 在 祁连山 上 不 同 海拔 种 群 的 遗传 多 样 性 水 平 的 研究 结论 ( 郭 
敏 等 ，2012)， 推 测 不 同 的 地 理 位 置 、 生 存 环境 和 海拔 高 度 都 体现 山 生 柳 不 显著 的 分 子 水 平 的 差 
异性 ， 但 是 差异 性 不 会 太 大 。 


3.2 山 生 柳 的 居 群 遗传 结构 


AMOVA 分 子 变 异 结果 表明 ， 山 生 柳 居 群 内 的 遗传 变异 占 总 变异 的 93.70% (CP<0.01)， 居 
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群 间 的 遗传 变异 仅 为 6.30%， 揭 示 了 山 生 柳 居 群 内 的 遗传 变异 比 居 群 间 的 遗传 变异 要 大 ， 居 群 
内 的 遗传 变异 是 山 生 柳 遗传 变异 的 主要 来 源 。 推 测 的 原因 : (1〉35 个 单 倍 型 中 仅 有 H01 12 
个 居 群 都 拥有 的 共享 单 倍 型 ， 共 有 208 个 个 体 拥有 ， 所 占 比例 为 59.26%， 由 于 H01 在 居 群 中 是 
无 规律 分 布 , 因此 山 生 柳 居 群 间 的 遗传 差异 并 不 明显 ;(2) 35 个 单 倍 型 中 有 22 个 单 倍 型 是 各 个 
居 群 所 独 有 ， 甚 至 某 一 居 群 会 同时 固定 了 多 个 特有 单 倍 型 ， 这 就 说 明山 生 柳 的 某 些 居 群 内 部 会 
有 可 辨别 的 遗传 变异 .特有 单 倍 型 比例 较 高 的 物种 往往 不 能 揭示 出 明显 的 谱系 地 理 格局 (Geng et 
al., 2018)。 繁 育 系 统 和 遗传 变异 分 布 之 间 存 在 着 密切 的 联系 CHogbin & Peakall, 1999)， 即 异 交 
物种 的 遗传 变异 大 多 分 布 在 群体 之 内 ， 和 群体 间 的 遗传 变异 通常 占 27% 以 下 《〈 赵 景泰 等 ，2008 )， 
与 本 研究 结果 相符 。 
植物 居 群 的 遗传 结构 受 群体 进化 历史 、 基 因 流 等 影响 ， 群 体 进化 历史 一 般 通 过 遗传 分 化 指 
数 CFsr) 来 反映 ， 其 大 小 可 在 一 定 程度 上 揭示 种 群 间 基因 流 和 和 遗传 漂 变 的 程度 程度 〈 杨 艳 婷 等 ， 
2018)。 一 般 来 讲 ， 当 Fsre [0，0.05)， 群 体 间 没有 遗传 分 化 ， 当 Fsre [0.05，0.15)， 群 体 分 
化 程度 中 等 ， 当 Fsre [0.15，0.25)， 群 体 间 高 度 分 化 ;Fsre [0.25，1]， 群 体 间 分 化 程度 非常 
5 (Wright, 1950; Wright，1970)。 山 生 柳 自然 居 群 的 遗传 分 化 指数 (Fsr) 为 0.06298， 表 明 
居 群 分 化 程度 为 中 等 。 山 生 柳 各 居 群 间 的 遗传 分 化 指数 (Fsr) 有 19 个 负 值 ， 占 28.79% CEW 
情况 下 Fsre [0，1])， 表 明 少 部 分 居 群 间 没 有 分 化 ， 但 总 体 上 还 是 反映 出 是 大 部 分 居 群 间 出 现 
非常 普遍 的 基因 交流 CWright, 1950; Wright, 1970; Wright, 1982)， 其 中 QM 和 NQ 之 间 的 Fsr 值 
最 大 为 0.2550。 山 生 柳 一 般 为 有 性 生殖 ， 花 粉 量 大 随 风 球 散 ， 但 在 严酷 的 生存 环境 下 会 进行 无 
性 繁殖 ， 本 研究 所 选 山 生 柳 材 料 地 理 分 布 较 广 、 地 理 环 境 高 度 复杂 ， 导 致 居 群 间 遗 传 变异 的 差 
异性 较 大 ， 有 的 居 群 几乎 没有 遗传 分 化 ‘QM, DT)， 有 的 居 群 分 化 程度 非常 高 CQL), 但 总 体 还 
是 表现 为 分 化 程度 中 等 偏 低 。 
根据 基因 流 公 式 ， 算 得 Nm=7.438， 当 Nm=7.438 > 1 ， 属 于 较 高 水 平 (Wright，1982)， 
c 表明 不 同 的 居 群 之 间 存 在 频繁 的 基因 流动 ， 增 大 群体 间 遗 传 差 异 的 表现 一 般 为 自然 选择 以 及 遗 
o 传 漂 变 ， 而 基因 流 可 以 “弱化 "群体 之 间 不 易 察 觉 的 遗传 差异 使 居 群 间 趋 向 统一 。 
对 山 生 柳 单 倍 型 变异 的 地 理 结构 检验 结果 为 Gst 值 为 0.068 和 Nsr [R73 0.075CNsr > Gs P > 
0.05)， 表 明山 生 柳 不 同居 群 间 出 现 了 大 量 的 亲缘 关系 靠近 的 单 倍 型 ， 可 以 证 明山 生 柳 不 存在 明 
显 的 谱系 地 理学 的 关系 。 基 于 遗传 距离 和 单 倍 型 的 UPGMA 聚 类 分 析 ， 得 到 山 生 柳 的 12 个 居 群 
各 分 为 四 组 且 与 12 个 居 群 地 理 分 布 没有 明显 相关 性 , 进一步 验证 了 山 生 柳 居 群 不 具有 谱系 地 理 
结构 。 


3.3 山 生 柳 的 种 群 扩张 

基因 流 分析 结 果 证 实 山 生 柳 的 12 个 居 群 存在 频繁 的 基因 交流 , 山 生 柳 居 群 间 不 存在 明显 的 
因 流 限制 。 通 过 对 所 有 地 理 居 群 单 倍 型 进行 失 配 分 析 ， 所 有 的 山 生 柳 单 倍 型 呈现 只 有 一 个 明 
顶峰 的 单 峰 型 曲线 ， 表 示 不 能 拒绝 居 群 扩张 的 假说 。 中 性 检验 的 结果 表明 ，Tajima’s D 
-2.28670, P < 0.01) 和 Fu's Fs (-5.29805, P < 0.020. 均 为 显著 的 负 值 ， 这 与 失 配 分 析 的 单 峰 型 
1 线 结果 相 一 致 。TCS 网 络 图 显示 HO 单 倍 型 为 祖先 单 倍 型 ， 由 于 H01 为 种 群 的 共享 单 倍 型 ， 


ng 
=$ 


Ko Bu 


故 不 能 确定 山 生 柳 种 群 的 扩散 中 心 。TCS 网 络 图 呈 星 状 辐射 分 布 ， 且 新 衍生 出 的 单 倍 型 与 原始 
单 倍 型 组 成 一 个 “ 星 状 ”辐射 结构 ， 说 明 该 群体 曾经 在 遇 到 瓶颈 后 迅速 扩 增 (Arbogast, 2000). 
通过 单 倍 型 多 样 性 与 核 苷 酸 多 样 性 分 析 得 到 单 倍 型 多 样 性 高 CHd > 0.50 MAEZ FEHER Pi 
< 0.005)， 预 示 着 群体 曾经 历 过 瓶颈 效应 后 ， 伴 随 了 迅速 的 种 群 扩张 与 变异 的 积累 。 综 上 所 述 ， 
推测 山 生 柳 在 其 分 布 范围 或 者 其 个 体 数 近 期 经 历 过 扩张 。 
对 青藏 高 原 地 区 相同 生境 的 高 山 绣 线 菊 (Khan et al., 2014)、 金 露 梅 (Ma et al., 2014). ib 
ee 
组 曾经 历 过 快速 地 种 群 扩张 事件 。 这 与 山 生 柳 的 种 群 扩张 结果 一 致 ， 推 测 物种 生境 和 繁育 方式 
的 相似 性 会 出 现 物种 为 了 营养 .繁殖 、 保 护 、 扩 散 以 及 对 环境 的 适应 出 现 协同 进化 的 现象 (Morgan 
et al., 2005; Thompson & Burdon, 1992)。 有 研究 表明 山 生 柳 种 群 可 能 以 有 性 和 无 性 两 种 繁殖 方式 
进行 更 新 以 巩固 其 优势 地 位 〈 李 妆 等 ，2002)， 年 龄 结构 分 布 图 呈正 三 角形 ， 种 群 动态 表现 为 增 
长 型 〈 李 毅 等 ，2002)， 这 与 本 研究 中 居 群 近期 出 现 过 扩张 相 吻合 。 花 粉 和 种 子 的 扩散 是 自然 植 
物种 群 最 主要 的 基因 流 方 式 〈 李 海 生 和 陈 桂 珠 ，2004)， 故 山 生 柳生 殖 方式 的 多 样 化 会 导致 其 基 
因 交流 频繁 ， 预 期 未 来 山 生 柳 居 群 会 处 于 进一步 扩张 状态 。 

山 生 柳 频 繁 的 基因 交流 和 经 历 过 近期 扩张 ， 猜 测 此 结果 可 能 的 原因 为 : (1) cpDNA 是 母系 
遗传 ， 估 算 的 基因 流 大 小 与 花粉 传播 无 关 ， 与 其 多 样 化 的 繁育 方式 有 关 ; (D 山 生 柳 一 般 以 有 
性 生殖 为 主 ， 其 生存 环境 是 长 期 积 雪 (低温 )、 干 燥 风 大 等 气象 特点 ， 使 其 种 子 可 以 随处 传播 ; 
(3) 由 于 种 子 的 传播 距离 远 ， 故 青藏 高 原 地 区 和 横断 山脉 错综复杂 地 形 也 只 能 在 很 小 的 程度 上 
阻隔 山 生 柳 居 群 间 的 基因 交流 ; (AO 在 严酷 环境 压力 下 ， 山 生 柳 种 子 活 力 极 低 ， 幼 苗 补充 量 受 
限 ， 从 而 有 性 生殖 受 限 ， 此 时 主要 依托 无 性 繁殖 来 巩固 其 物种 数量 的 稳定 〈 陈 文 业 等 ，2008 )。 
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